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We consider the classical and quantum properties of the " Chirikov typical map" , proposed by Boris 
Chirikov in 1969. This map is obtained from the weU known Chirikov standard map by introducing 
a flnite number T of random phase shift angles. These angles induce a random behavior for small 
time scales (t < T) and a T-periodic iterated map which is relevant for larger time scales {t > T). 
We identify the classical chaos border kc ~ j—s/z ^ ^ for the kick parameter k and two regimes 
with diffusive behavior on short and long time scales. The quantum dynamics is characterized by 
the effect of Chirikov localization (or dynamical localization). We And that the localization length 
depends in a subtle way on the two classical diffusion constants in the two time-scale regime. 

PACS numbers: 05.45.Mt,05.45.Ac,72.15.Rn 
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I. INTRODUCTION 



The dynamical chaos in Hamiltonian systems often 
leads to a relatively rapid phase mixing and a relatively 
slow diffusive spreading of particle density in an action 
space A well known example of such a behavior is 
given by the Chirikov standard map [2,, ,3,] . This simple 
area-preserving map appears in a description of dynamics 
of various physical systems showing also a generic behav- 
ior of chaotic Hamilton systems [4] . The quantum version 
of this map, known as the quantum Chirikov standard 
map or kicked rotator, shows a phenomenon of quan- 
tum localization of dynamical chaos which we will call 
the Chirikov localization, it is also known as the dynam- 
ical localization. This phenomenon had been first seen 
in numerical simulations [5] while the dependence of the 
localization length on the classical diffusion rate and the 
Planck constant was estabhshed in [6, 7, 8]. It was shown 
in that this localization is analogous to the Ander- 
son localization of quantum waves in a random potential 
(see [lo| for more Refs. and details). In this respect 
the Chirikov localization can be viewed as a dynamical 
version of the Anderson localization: in dynamical sys- 
tems diffusion appears due to dynamical chaos while in 
a random potential diffusion appears due to disorder but 
in both cases the quantum interference leads to localiza- 
tion of this diffusion. The quantum Chirikov standard 
map has been built up experimentally with cold atoms 
in kicked optical lattices . In one dimension all states 
remain localized. In systems with higher dimension (e.g. 
(i = 3) a transition from localized to delocalized behavior 
takes place as it has been demonstrated in recent exper- 
iments with cold atoms in kicked optical lattices [l2| . 

While the Chirikov standard map finds various appli- 
cations it still corresponds to a regime of kicked systems 
which are not necessarily able to describe a continue flow 
behavior in time. To describe the properties of such a 
chaotic flow Chirikov introduced in 1969 [2,] a typical map 
which we will call the Chirikov typical map. It is obtained 
by T iterations of the Chirikov standard map with ran- 



dom phases which are repeated after T iterations. For 
small kick amplitude this model describes a continuous 
flow in time with the Kolmogorov- Sinai entropy being 
independent of T. In this way this model is well suited 
for description of chaotic continuous flow systems, e.g. 
dynamics of a particle in random rnagnetic fields [l^] or 
ray dynamics in rough billiards [l^^Illl- Till present only 
certain properties of the classical typical map have been 
considered in [H, @ . 

In this work we study in detail the classical diffusion 
and quantum localization in the Chirikov typical map. 
Our results confirm the estimates presented by Chirikov 
p!] for the classical diffusion I^o and instability. For the 
quantum model we find that the localization length i is 
determined by the classical diffusion rate per a period of 
the map £ ~ DqT/H^ in agreement with the theory 

The paper is constructed in the following way: Section 
II gives the model description, dynamics properties and 
chaos border are described in Section III, the properties 
of classical diffusion are described in Section IV, the Lya- 
punov exponent and instability properties are considered 
in Section V, the quantum evolution is analyzed in Sec- 
tion VI, the Chirikov localization is studied in Section 
VII and the discussion is presented in Section VIII. 



II. MODEL 

The physical model describing the Chirikov typical 
map can be obtained in a following way. We consider 
a kicked rotator with a kick force rotating in time (see 
Fig. IB): 



F(t) = k 



- cos a{t) 
sin a{t) 



(1) 



where /c is a parameter characterizing the amplitude of 
the kick force. We measure the time in units of the ele- 
mentary kick period and we assume that the angle a{t) 
is a T-periodic function a{t + T) = a{T) with T being 
integer (i. e. an integer multiple of the elementary kick 
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period). The time dependent Hamiltonian associated to 
this type of kicked rotator reads: 



m) 



■kcos[0 ^a{t)] S{t-n). (2) 




FIG. 1: Kicked rotator where the kick force rotates in time 
with a kick force angle a{t) that is a periodic function in time. 

To study the classical dynamics of this Hamiltonian 
we consider the values of p and slightly before the kick 
times: 



Pt 



lim p(t 

s\0^ 



lim 0(t 

s\0 



where t is the integer time variable. The time evolution 
is governed by the Chirikov typical map defined as: 



Pt+i =Pt^k sm{Ot + at) 



(3) 



with at = a{t) at integer times. Since at = at mod t 
there are T independent different phase shifts at for t G 
{0, T — 1}. Following the approach of Chirikov [If 
we assume that these T phase shifts are independent and 
uniformly distributed in the interval [0,27r[. 

For the quantum dynamics we consider the quantum 
state slightly before the kick times: 



lim \ij{t - 

£\0 



■e)) 



whose time evolution is governed by the quantum map: 



IV^t+i) = exp 



.p' , 



k 

-i— cos{6 



(4) 

with p = —ihd/dO and the wave function ilj{0 + 27r) = 
iIj{6) = {6\iIj). In the following we refer to the map (|3]), 
(HI as the classical or quantum version of the Chirkov 
typical map as it was introduced by Chirikov [2]. The 
Chirikov standard map [3] corresponds to a particular 
choice of random phases all being equal to a constant. 



III. CLASSICAL DYNAMICS AND CHAOS 
BORDER 

It is well known that the Chirikov standard map ex- 
hibits a transition to global (diffusive) chaos at /c > /Cc ~ 



0.9716 [H. For the Chirikov typical map ^ it is pos- 
sible to observe the transition to global chaos at values 
/Cc <C 1 provided that T 1. In order to determine kc 
quantitatively we apply the Chirikov criterion of overlap- 
pmg resonances [16] (see also [i,[3| for more details). For 
this we develop the T-periodic kick potential in ([2|) in the 
Fourier series: 



^i{6 —mujt) 



(5) 



with: 



rj, •) Jm — J, / ^ ^ 



(6) 



The resonances correspond to 6{t) = muot with integer m 
and positions in p-space: Pres = = ^iim/T. There- 
fore the distance between two neighboring resonances is 
Ap = uo = 27r/T and the separatrix width of a resonance 
is - 4(|/^|2)V4 = 4 (^2/^)1/4 ^^^^^ 

age is done with respect to a^. The transition to global 
chaos with diffusion in p takes place when the resonances 
overlap, i. e. if their width is larger than their distance 
which corresponds to the chaos border 



k>kc = 7rV(4T^/^) 



(7) 



with A:c < 1 for T > 1. 



k = 0.05 



k - 0.078 



A; = 0.1 




FIG. 2: Classical Poincare sections of 50 trajectories of length 
tmax = 10000 of the Chirikov typical map for one particular 
realization of the random phase shifts with T = 10 and three 
values: k = 0.05, k = 0.078, k = 0.1. The positions (Ot.pt) G 
[0, 27r[x [0, 27r[ are shown after applications of the full map, 
i. e. T-times iterated typical map with: t mod T = 0. The 
value of the middle figure corresponds to the critical value 
kc ~ 0.078 of dZl) for the transition to global diffusive chaos. 

It is well known [l|, [H, 0] that the Chirikov criterion of 
overlapping resonances gives for the Chirikov standard 
map a numerical value 7r^/4 2.47 which is larger than 
the real critical value ~ 0.9716 that is due to a combina- 
tion of various reasons (effect of second order resonances, 
finite width of the chaotic layer etc.). However, for the 
Chirikov typical map these effects appear to be less im- 
portant, probably due of the random phases in the ampli- 
tudes fm and phase shifts of resonance positions. Thus, 
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the expression ([7j) for the chaos border works quite weh 
including the numerical prefactor, even though the ex- 
act value kc depends also on a particular realization of 
random phases. In Fig. [2] we show the Poincare sec- 
tions of the Chirikov typical map for a particular ran- 
dom phase realization for T = 10 and three different 
values of /c G {0.05,0.078,0.1} where the critical value 
at T = 10 is kc = (f )2 10-^/^ ^ q Q^g Fig. [2] confirms 
quite well that the transition to global chaos happens at 
that value. Actually choosing one particular initial con- 
dition ^0 = 0-8 • 27r and po = 0.25 • 27r one clearly sees 
that at k = 0.05 only one invariant curve at p po is 
filled and that at /c = 0.1 nearly the full elementary cell 
is filled diffusively. At the critical value k = 0.078 only 
a part of the region p > po is diffusively filled during a 
given number of map iterations. 



IV. CLASSICAL DIFFUSION 

For k > kc the classical dynamics becomes diffusive in 
Pt and for k ^ kc we can easily evaluate the diffusion 
constant assuming that the angles 6t are completely ran- 
dom and uncorrelated. In order to discuss this in more 
detail we iterate the classical map (|3]) up to times t: 

t-1 

Pt = Po ^ k^ sin(6>i. + a^.) (8) 

1^=0 

and 

t t-1 
Ot = Oo^'^p^ = Oo^pot^k'^it-iy) sin((9^+a^) . (9) 

/i=l u=0 

Using the assumption of random and uncorrelated angles 
we find that the quantities 6pt = Pt — Po and 66t = Of — 
Oo — Pot are random gaussian variables (for t ^ 1) with 
average and variance: 

(Spt) = {S0t) = , (10) 

(11) 

implying a diffusive behavior in p-space with diffusion 
constant I^o = (Ap)^/At = k'^/2 which is valid for k > 
kc. For k > kc but close to kc we expect the dynamics 
also to be diffusive but with a reduced diffusion constant 
D < Do due to correlations of Ot for different times t and 
for k < kc we have = 0. However, we insist that this 
behavior is expected for long times scales, in particular 
t ^ T with T being the period of the random phases 
o^u = OLT+v For small times t < T there is always, for 
arbitrary values of k (including the case k <kc)^ ^ simple 
short time diffusion with the diffusion constant I^o and 
in this regime Eq. ([TT]) is actually exact (if the average 
is understood as the average with respect to a^). 



It is interesting to note that Eq. ([TT]) provides an ad- 
ditional way to derive the chaos border kc. Actually, 
we expect the long term dynamics to be diffusive if at 
the end of the first period t = T the short time diffu- 
sion allows to cross at least one resonance in p space : 
8pT ^ kVr > Ap = 27r/T or if the dynamics becomes 
ergodic in space : SOt - feT^/^ > 27r. Both condi- 
tions provide the same chaos border k > kc ^ 
exactly confirming the finding of the previous section by 
the Chirikov criterion of overlapping resonances. 

We note that the behavior 50t ~ k t^l^ is a direct conse- 
quence that bdi is a sum (integral) over bp^ for /j. < t and 
that Pj^ itself is submitted to a diffusive dynamics with 
Sp^ ^ ji^^'^. The same type of phase fluctuations also 
happens in other situations, notably in quantum dynam- 
ics where the quantum phase is submitted to some kind 
of noise with diffusion in energy or in frequency space 
(with diffusion constant I^o) implying a dephasing time 
t(j) ^ Dq . For example, a similar situation appears 
for dephasing time in disordered conductors [13] where 
the diffusive energy fluctuations for one-particle states 
are caused by electron-electron interactions and where 
the same type of parametric dependence (in terms of the 
energy diffusion constant) is known to hold. Another 
example is the adiabatic destruction of Anderson local- 
ization discussed in [18] where a small noise in the hop- 
ping matrix elements of the one-dimensional Anderson 
model leads to a destruction of localization and diffusion 
in lattice space because the quantum phase coherence is 
limited by the same kind of mechanism. 

We now turn to the discussion of our numerical study 
of the classical diffusion in the Chirikov typical map. We 
have numerically determined the classical diffusion con- 
stant for values of 10 < T < 1000 and T-^/^ < k < 1. 
For this we have simulated the classical map up to times 
t < 100 T and calculated the long time diffusion con- 
stant D as the slope from the linear fit of the variance 
(Sp^) for lOT < t < 100 T (in order to exclude artifi- 
cial effects due to the obvious short time diffusion with 
1^0 = k'^/2). According to Fig. [3] the ratio D/Dq can be 
quite well expressed as a scaling function of the quantity : 
kT^/^ - k/kc: 

D ^ Dof (/cT^/') (12) 

where the scaling function f{x) can be approximated by 
the fit : 



/,.)=»p^l-/l + l|? + H|l]^ ,13) 

This fit has been obtained from a plot of \ii{D/Do) ver- 
sus y = where the numerical data gives a linear 
behavior \n{D/Do) ^ -Aiy for < 1 and \n{D/Do) ^ 
— 742^+const. for y 1 with different slopes Ai and 
A2 which can be fitted by the ansatz \ii{D/Do) = 1 — 
^/l -\- Ci y ^ C2 y'^ . The problem to obtain an analytical 
theory for this scaling function is highly non-trivial and 
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FIG. 3: (Color online) The ratio D / Dq of the classical diffu- 
sion constant D of the Chirikov typical map over the theo- 
retical diffusion constant Dq — assuming perfectly un- 
correlated phases ^t, as a function of the scaling parameter 
X = /cT^/^ for 10 < T < 1000 and T'^^^ < k < 1. The 
classical diffusion constant D has been obtained from a lin- 
ear fit of the average variance < 6pt > in the time interval 
lOT < t < 100 T. The average variance has been calculated 
for 400 different realizations of the random phase shifts and 
25 different random initial conditions (^o,Po) ^ [0, 27r[x [0, 27r[ 
for each realization. The vertical full (blue) line represents 
the classical chaos border at kcT^^^ — 7r^/4 ^ 2.47 (compare 
Fig. [2]) and the grey rectangle at the left shadows the non- 
diffusive regime where the numerical procedure (incorrectly) 
yields small positive values of D (see text). The full curved 
(black) line represents the scaling function D / Dq — f(x) as 
given by Eq. (|13|) . The inset shows DT^ as a function of the 
scaling parameter x = kT^^'^ and the full (red) line represents 
the scaling function x^/(x)/2. The full vertical (blue) line 
again represents the classical chaos border. 



subject to future research. However, we note that the nu- 
merical scaling function (fT3|) shows the correct behavior 
in the limits x ^ oo and x ^ 0. For x = k T^/^ > IQQ we 
have f{x)^l implying D ^ Dq for the strongly diffusive 
regime where phase correlations in Of can be neglected. 
For intermediate values kcT^^'^ ^ 2.47 < x < 100 there 
is a two scale diffusive regime with a short time diffusion 
constant I^o for t < T and a reduced longer time diffusion 
constant D < Do for t > T. 

For k < kc {x < 2.47) the long time diffusion con- 
stant is expected to be zero but the numerical fit proce- 
dure still results in small positive values (note that the 
scaling function vanishes very quickly in a non-analytical 
way f{x) ^ exp(-4.8/x) as X 0) simply because here 
the chosen fit interval lOT < t < 100 T is too small. 
In order to numerically identify the absence of diffusion 
it would be necessary to consider much longer iteration 
times. We also note that the variance {Sp^) shows a quite 
oscillatory behavior for k < kc and t > T indicating the 
non-diffusive character of the dynamics despite the small 
positive slope which is obtained from a numerical fit (in 



a limited time interval). Therefore, we have shadowed in 
Fig. [3] the regime x < 2.47 by a grey rectangle in order 
to clarify that this regime is non- diffusive. 

The important conclusion of Fig. [3] is that it clearly 
confirms the transition to chaotic diffusive dynamics at 
values k > kc ^ 2^-3/2 ^^^^ addition there is 
even an approximate scaling behavior in the parameter 
X = kT^/^ - k/kc. 

V. LYAPUNOV EXPONENT AND ERGODIC 
DEPHASING TIME SCALE 

We consider the chaotic regime kc < k < 1 where the 
diffusion rate is quite slow and where we expect that the 
Lyapunov exponent is much smaller than unity implying 
that the exponential instability of the trajectories devel- 
ops only after several iterations of the classical map. In 
order to study this in more detail, we rewrite the map 
as : 

Pt+i = + MOt) , 6>t+i = Ot +pt+i (14) 

where for the Chirikov typical map we have ft{0) = 
kcos{0 + at) but in this section we would like to al- 
low for more general periodic kick functions ft{0) with 
vanishing average (over 0). In order to determine the 
Lyapunov exponent of this map we need to consider two 
trajectories iOt^Pt) and iOt^Pt) both being solutions of 
(p!4|) with very close initial conditions at t = 0. The dif- 
ferences Apt = Pt — Pt and AOf = Ot — Ot are iterated by 
the following linear map: 

Apt^i ^ Apt + fi{0t)A0t , AOt^i = AOt + Ap,+i 

(15) 

as long as |A^^| <C 1. In a similar way as with the Eqs. 
dH]) and ([9]) in the last section, we may iterate the linear 
map up to times t : 

t-i 

Apt = Apo + Y^fl{e^,)Ae, , (16) 

t 

AOt = A^o + ^Ap^ 

t-i 

= AOo^tApo^Y.(t-u)fUO,)AO,. (17) 

We now assume that in the chaotic regime the phases 0^ 
are random and uncorrelated, and that the average of the 
squared phase difference behaves as : 

{Ae^) = e{t)A9l (18) 

with a smooth function S{t) we want to determine. From 
Eq. ([T7|) we obtain in the continuum limit the following 
integral equation for the function Q{t) : 

G{t) = 1^2At^A^t'^^K [ {t-vfe{v)dv (19) 
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where k, = {fl{0)'^)e and A = Apq/AOq. For the map 
(|3]) we have tz = This integral equation imphes the 

differential equation 

&'\t) =2f^e{t) (20) 
with the general solution : 

3 

e{t) = ^Cje^^^' (21) 

where the constants Cj are determined by the initial 
conditions at t = and 2Xj are the three solutions of 
{2Xj)^ = 2i<i : 

Ai = (2/^)1/^2 , A2 = e'"'/'Ai , A3 = . (22) 

Since only Ai has a positive real part, we have in the long 
time limit : 

Q(t)^Cie^^'' (23) 

and therefore by Eq. ([18]) Ai represents approximately 
the Lyapunov exponent of the map (p!4|) . 

We note that this calculation of the Lyapunov expo- 
nent is not exact, essentially because we evaluate the di- 
rect average (AO^) instead of exp((ln(A^^))). A proper 
and more careful evaluation of the Lyapunov exponent 
for this kind of maps (in the regime <C 1, assuming 
chaotic behavior with uncorrelated phases) has been done 
by Rechester et al. [l3| in the context of a motion along 
a stochastic magnetic field. Their result reads in our no- 
tations : 

A = — ^ 0-36 ^ 0.29^^^ (24) 

while from Eq. ([22]) we have Ai ^ 0.63 a^^/^ with the 
same parametric dependence but with a different numer- 
ical prefactor. The dependence ([24]) has been confirmed 
in numerical simulations [6] for the model ([3]) and we do 
not perform numerical simulations for A here. We note 
that in the map ([3]) the Lyapunov exponent A gives the 
Kolmogorov- Sinai entropy [1]. 

The inverse of the Lyapunov exponent defines the 
Lyapunov time scale ^Lyap ^ k,~^^^ which is the time 
necessary to develop the exponential instability of the 
chaotic motion. According to Eq. ([TT]) we also have : 
{60^) = Do t^/3 with I^o = {fu{0)'^)e implying an ergodic 
dephasing time ^ Dq being the time necessary for 
a complete dephasing where there is no correlation of 
the actual phase Ot with respect to the ballistic phase 
^baii. = -\- Pot. For the Chirikov typical map with 
fu{0) = k s'm{0 + au) the averages of fl{0)'^ and fviO^ 
are identical, implying n = Dq = /c^/2, and these two 
time scales coincide : ^Lyap = ^ k~'^^^. Furthermore, 
the condition for global chaos, k > kc ^ 2^-3/2^ reads 
^Lyap = < T implying that the exponential instability 



and complete dephasing must happen before the period 
T. 

We mention that for other type of maps, in particular 
if fu{0) contains higher harmonics such as sin(M^) we 
may have : ~ M'^Dq, and therefore parametrically 
different times scales ^Lyap ~ t^. Examples of these 
type of maps have been studied in Ref. [l3| and also 
in Refs. [3, llH in the context of angular momentum 
diffusion and localization in rough billiards. 

For the Chirikov typical map there is a further time 
scale ^Res which is the time necessary to cross one res- 
onance of width Ap = 27r/T by the diffusive motion 
{5pl) = kH/2: 

(25) 

In the chaotic regime k > kc this time scale is parametri- 
cally smaller than the dephasing time and the Lyapunov 
time scale. 



VI. QUANTUM EVOLUTION 

We now study the quantum evolution which is de- 
scribed by the quantum map ([4]) (see section |II]) . Typ- 
ically in the literature studying the quantum version of 
Chirikov standard map the value of h is absorbed in a 
modification of the elementary kick period and the kick 
parameter k. Here we prefer to keep h as an indepen- 
dent parameter (also for the numerical simulations) and 
to keep the notation T of the (integer) period of the time 
dependent kick-angle a{t). In this way we clearly iden- 
tify two independent classical parameters k and T, one 
quantum parameter h and a numerical parameter N ^ 1 
being the finite dimension of the Hilbert space for the 
numerical quantum simulations. For the physical under- 
standing and discussion it is quite useful to well separate 
the different roles of these parameters and we further- 
more avoid the need to translate between "classical" and 
"quantum" versions of the kick strength k. 

We choose the Hilbert space dimension to be a power 
of 2 : N = 2^ allowing an efficient use of the discrete 
fast Fourier transform (FFT) in order to switch between 
momentum and position (phase) representation. A state 
\ip) is represented by a complex vector with N elements 
'^(j)? j = 0, . . . , — 1 where pj = hj are the discrete 
eigenvalues of the momentum operator p. In order to 
apply the map ([4]) to the state |?/;) we first use an inverse 
FFT to transform to the phase representation in which 
the operator 6 is diagonal with eigenvalues 6j = 27rj/N, 
j = 0,...,A/' — 1, then we apply the first unitary matrix 
factor in ([4]) which is diagonal in this base, we apply an 
FFT to go back to the momentum representation and 
finally the second unitary matrix factor, diagonal in the 
momentum representation, is applied. This procedure 
can be done with 0{N \og2{N)) operations (for the FFT 
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and its inverse) plus 0{N) operations for the application 
of the diagonal unitary operators. 

We have also to choose a numerical value of h and this 
depends on which kind of regime (semiclassical regime or 
strong quantum regime) we want to investigate. We note 
that due to the dimensional cut off we obtain a quantum 
periodic boundary condition : V^(0) = '^(A^), i. e. the 
quantum dynamics provides in p representation always 
a momentum period of hN. Furthermore the classical 
map is also periodic in momentum with period 27r and 
it is possible to restrict the classical momentum to one 
elementary cell where the momentum is taken modulo 27r. 
Therefore one plausible choice of h amounts to choose the 
quantum period hN to be equal to the classical period 
27r, i. e.: h = 2t:/N. In this case the quantum state covers 
exactly one elementary cell in phase space. Since ^ 
for TV ^ oo, we refer to this choice as the semiclassical 
value of h. 

In this section, we present some numerical results of the 
quantum dynamics using the semiclassical value of h. In 
Fig.HJ we show the Husimi functions for the case T = 10, 
k = 0.1 with an initial state being a minimal gaussian 
wave packet centered at 6^0 = 0.8 • 27r and po = 0.25 • 27r 
with a variance (in ^-representation) of Ap = 27r/Vl27V. 
This position is well inside a classically chaotic region 
(see Fig. [2]). The Hilbert space dimensions N are 2^^ and 
2^^ and the iteration times are t G {0, 20, 60, 100, 150}. 
We see that at these time scales the motion extends to 
two classical resonances with two stable and quite large 
islands. For N = 2^'^ the classical phase space structure 
is quite visible but the finite "resolution" in phase space 
due to quantum effects is quite strong while for N = 2^^ 
the Husimi function allows to resolve much better smaller 
details of the classical motion. We note that the Husimi 
function is obtained by smoothing of the Wigner function 
over a phase space cell of h size (see e.g. more detailed 
definitions and Refs. in [l^). 

We note that for T = 10 the value of the kick strength 

= 0.1 is only slightly above the chaos border 0.078 with 
global diffusion but there are still large stable islands that 
occupy a significant fraction of the phase space. At this 
value the numerically computed diffusion constant is D = 
OADo ^ 1/2000 (see section [Hi] and Fig. In Fig. [5]we 
compare the Husimi functions for various Hilbert space 
dimensions at time t = 20000. This time is sufficiently 
long so that a diffusive spreading in p-space gives 6pt = 
\fDi ^ 3.162 thus roughly covering one elementary cell 
(but not absolutely uniformly). 

For the smallest value of = 2^ we observe a very 
strong influence of quantum effects with a Husimi func- 
tion extended to half an elementary cell which is signif- 
icantly stronger localized (in the momentum direction) 
than the diffusive classical spreading would suggest. Fur- 
thermore we cannot identify any classical phase space 
structure. This is quite normal due to the very limited 
resolution of \/N = 16 "quantum-pixels" in both direc- 
tions of and p. For N = 2^^ the quantum effects are 
still strong but the Husimi function already extends to 




FIG. 4: (Color online) Quantum evolution of the Chirikov 
typical map at k = 0.1 and T = 10 in the semiclassical limit 
with h — 2'k/N and with the Hilbert space dimension N — 2^'^ 
(left column) and N — 2^^ (right column) at times t — (first 
row), t = 20 (second row), t = 60 (third row), t — 100 (fourth 
row) and t — 150 (fifth row). Shown are Husimi functions 
with maximum values at red (grey), intermediate values at 
green (light grey) and minimum values at blue (black) at the 
lower half of one elementary cell G [0,27r[, p G [0,7r[. The 
initial condition is a coherent gaussian state centered at Oq = 
0.8 •27r and po = 0.25 •27r with a variance (in ^-representation) 
of Ap = 2n / y/V2N . The resolution corresponds to ^/N (64 or 
256) squares in one line. The realization of the random phase 
shifts is identical to that of Fig. [2] Note that at the considered 
value k — 0.1 the global classical dynamic is diffusive but 
requires iteration times of t ~ 10000 to fill one elementary 
cell (see Fig. [5]). 



85% of the elementary cell and we can identify first very 
slight traces of at least three large stable islands. For 
TV = 2^^ the Husimi function fills the elementary cell 
as suggested by the classical spreading but the distribu- 
tion is less uniform than for larger values of N or for the 
classical case. We can also quite well identify the large 
scale structure of the phase space with the main stable 
islands associated to each resonance. However, the fine 
structure of phase space is not visible due to quantum 
effects. For N = 2^^ and even more for TV = 2^^, the res- 
olution of the Husimi function increases and approaches 
the classical distribution which is also shown in Fig. [5] 
for comparison. For N = 2^^, we even see first traces 
of small secondary islands that are not associated to the 
main resonances. We note that the classical distribution 
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showing that the spearing in phase reaches 27r approxi- 
mately at t = 100. 

In summary the quantum simulation of the Chirikov 
typical map using the semiclassical value h = 27t/N re- 
produces quite well the classical phase space structure for 
sufficiently large N while for smaller values of the Hilbert 
space dimension {N < 2^) the nature of motion remains 
strongly quantum and the diffusive spreading over the 
cell is stopped by quantum localization. 



VII. CHIRIKOV LOCALIZATION 

It is well established 0, Q that, in general, quan- 
tum maps on one-dimensional lattices, whose classical 
counterpart is diffusive with diffusion constant Dc\^ show 
dynamical exponential localization of the eigenstates of 
the unitary map operator with the localization length 
io = Dci./Ji^ measured in number of lattice sites (that 
corresponds to the quantum number n associated to the 
momentum hy p = hn). Here I^ci. is the diffusion rate 
in action per period of the map. This expression for 
is valid for the unitary symmetry class which applies to 
the Chirikov typical map which is not symmetric with 
respect to the transformation 6^—6. We have further- 
more to take into account that the map (|4|) depends on 
time due to the random phases and in order to deter- 
mine the localization length we have to use the diffusion 
constant for the full T-times iterated map (which does 
not depend on time) : I^ci. = DqT assuming we are in 
the regime where D = I^o foi" k ^ kc. In this case we 
expect a localization length 



FIG. 5: (Color online) Quantum evolution of the Chirikov 
typical map as in Fig. |4] with the same coherent gaussian 
state as initial condition, same values /c = 0.1 and T = 10, 
h = 2ti /N but at the iteration time t — 20000 and at different 
values of AT = 2^ AT = 2^° (first row), AT = 2^^ AT = 2^^ 
(second row), N — 2^^ and classical simulation (third row). 
For the classical map 20000 trajectories have been iterated up 
to the same time t — 20000 with random initials conditions 
very close to the initial position at = 0.8 • 2n it 0.002 and 
po = 0.25 • 27r ± 0.002. Colors are as in Fig. H 



in Fig. [5] is not a full phase portrait (showing "all" iter- 
ation times) but it only contains the classical positions 
after t = 20000 iterations with random initial positions : 
l9o = 0.8 • 27r ± 0.002, po = 0.25 • 2it ± 0.002 close to the 
gaussian wave packet used as initial state in Figs. [H [5l 
Finally we note that a wave packet with initial size 
grows exponentially with time and spreads 
over the whole phase interval 2it after the Ehrenfest time 
scale 0, [lO, [m 



tE ^ ln(27r/V^)/A. 



(26) 



For the parameters of Fig. [H e.g. h = 27r/2^^, this gives 
tE ~ 103 that is in agreement with the numerical data 
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fi2 



(27) 



This localization length is obtained as the exponential 
localization length from the eigenvectors of the full (T- 
times iterated) unitary map operator. In numerical stud- 
ies of the Chirikov typical map it is very difficult and 
costly to access to these eigenvectors and we prefer to 
simply iterate the quantum map with an initial state lo- 
calized at one momentum value in ^-representation and 
to measure the exponential spreading of |'^(t)) at suffi- 
ciently large times (using time and ensemble average). 
This procedure is known Q to provide a localization 
length i artificially enhanced by a factor of two: i = 2£o. 

Let us note that the relation ([27j) assumes that the clas- 
sical dynamics is chaotic and is characterized by the diffu- 
sion 1^0 • It also assumes that the eigenvalues of the uni- 
tary evolution operator are homogeneously distributed 
on the unitary circle. For certain dynamical chaotic sys- 
tems the second condition can be violated giving rise to 
a multifractal spectrum and delocalized eigenstates (this 
is e.g. the case of the kicked Harper model [22|,[2^). The 
analytical derivation of (|271) using supersymmetry field 
theory assumes directly jl5| or indirectly that the 
above second condition is satisfied. We also assume that 
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this condition is satisfied due to randomness of at phases 
in the map (|3]). 

Before we discuss our numerical results for the local- 
ization length we would like to remind the phenomeno- 
logical argument coined in |6| which allows to determine 
the above expression relating localization length to the 
diffusion constant and which we will below refine in order 
to take into account the two-scale diffusion with different 
diffusion constants at short and long time scales. 

Suppose that the classical spreading in p-space is given 
by a known function : 



P2{t) = {5pf, 



(28) 



where P2(t) is a linear function for simple diffusion but 
it may be more general in the context of this argumenta- 
tion (but still below the ballistic behavior for t ^ oo). 
For the quantum dynamics we choose an initial state lo- 
calized at one momentum value. This state can be ex- 
panded using ^ eigenstates of the full map operator with 
a typical eigenphase spacing 27r/^. The iteration time t 
(with t being an integer multiple of T) corresponds to 
t/T applications of the full map operator. We expect 
the quantum dynamics to follow the classical spreading 
law (|28|) for short time scales such that we cannot re- 
solve individual eigenstates of the full map operator, i. e. 
for 27Tt/{T£) < 2iT. Therefore at the critical time scale 
t* = IT we expect the effect of quantum localization to 
set in and to saturate the classical spreading at the value 
^P\oc ^ ^^^^ finite localization length £ (mea- 

sured in integer units of momentum quantum numbers). 
This provides an implicit equation for t* : 



e = c 



P2(r) 



(29) 



where C is a numerical constant of order unity. Let us 
first consider the case of simple diffusion with constant 
1^0 for which we have P2it) = Dot- In this case we 
obtain : 



C 



do 



(30) 



with £o given by Eq. (|27j) . This result provides the nu- 
merically measured value £ — 2£o by the exponential 
spreading of |'0(t)) if we choose (7 = 2. Below we will 
also apply Eq. ([29|) to the case of two scale diffusion 
with different diffusion constants at short and long time 
scales providing a modified expression for the localization 
length. 

First we want to present our numerical results for the 
localization length of the quantum Chirikov typical map. 
Since the localization length scales as h~'^ we cannot 
use the semiclassical value h = 2it/N since in the limit 
N ^ oo the localization length would always be larger 
than N and in addition we would only cover one ele- 
mentary classical cell of phase space which is not very 
suitable to study the effects of diffusion and localization 
in momentum space. Therefore we choose a finite and 
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FIG. 6: (Color online) Illustration of the Chirikov localization 
for the quantum Chirikov typical map. Shown is the (aver- 
aged) absolute squared wave function in p-representation as 
a function of the integer quantum number n associated to 
p — hn. The initial state at t = is localized at n = 
and the quantum typical map has been applied up to times 
t — tmax with tmax choscu such that the initial classical dif- 
fusion is well saturated at tmax/4 due to quantum effects (see 
Fig. [7|. In order to determine numerically the localization 
length, first iV^nl^ has been time-averaged for the time interval 
tmax/4 < t < tmax and thcu a further ensemble average (over 
100 realizations of the random phase shifts) of \ipn\'^ (upper 
red/gray curve) or of Ind^^nl^) (lower blue/black curve) has 
been applied. The localization lengths are obtained as the 
(double) inverse slopes of a linear fit (with relative weights 
Wn iV^nl^) of ln(|'0n|^) vcrsus 71. Here the parameters are: 
N = 2^^ k = 0.5, T = 100, tmax = 546132, h = 27r/(17 + 7) ^ 
0.357 with the golden number 7 = ( ^5 - l)/2 ^ 0.618. The 
fits are given by the thick straight lines. The theoretical local- 
ization length is io = k^T/{2h^) ^ 98.3 while the numerical 
fits provide ~ 231.7 (upper curve) and i\n^ ~ 187.5 (lower 
curve) . 



fixed value h = 2it/{N + 7) where N is some fixed inte- 
ger in the range 1 <C TV <C iV and 7 = (\/5 - 1) /2 ^ 0.618 
is the golden number because we want to avoid artificial 
resonance effects between the classical momentum period 
27r and the quantum period HN (due to the finite dimen- 
sional Hilbert space). The ratio of these two periods, 
N/{N + 7) ^ 1, is roughly the number of elementary 
classical cells covered by the quantum simulation. For 
most simulations we have chosen = 17 and varied the 
classical parameters k and T but we also provide the data 
for a case where the values of k and T are fixed and N 
varies. 

For the numerical quantum simulation we choose the 
initial state |'0(O)) = |0) being perfectly localized in mo- 
mentum space and we apply the quantum Chirikov typ- 
ical map up to a sufficiently large time scale tmax chosen 
such that the initial diffusion is well saturated at tmax/4 



and we perform a time average of \ipn\ 



for 



the time interval tmax/4 < t < tmax- The resulting time 
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average is furthermore averaged with respect to different 
reaUzations of the T random phases aj, and here we con- 
sider two cases where we average either l^/^^P or Ind^/^^Pj- 
We then determine two numerical values £^ and £in ^ by a 
linear fit of ln(|?/^np) = 2n/^+const. with relative weight 
factors Wn ^ IV^nP in order to emphasize the initial ex- 
ponential decay and to avoid problems at large values of 
n where the finite numerical precision (~ 10~^^) or the 
finite value of N may create an artificial saturation of the 
exponential decay. This procedure is illustrated in Fig. [6] 
for the parameters k = 0.5, T = 100, h = 2i\ j (I7+7) and 
N = 2^^. We see that the two numerical values and 
^in^ are roughly 2^o as expected but may differ among 
themselves by a modest numerical factor. 
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FIG. 7: (Color online) The ensemble averaged variance of the 
expectation value of v? as a function of i for the same parame- 
ters (and the same realizations of the random phase shifts) as 
in Fig. [6] The upper horizontal (red) line corresponds to the 
saturation value ^^iff (obtained from a time average of < > 
in the interval tmax/4 < t < tmax) with -£diff ~ 146.8 being the 
localization due to saturation of diffusion. The short lower 
horizontal (green) line corresponds to with the theoretical 
localization length io ^ 98.3 (see also Fig. [G]). The inset 
shows the initial diffusive regime at shorter time scales. The 
straight (red/grey) line shows the value Dq/H^ ^ 0.983 asso- 
ciated to the initial diffusion: {in?) — {p^) ^ Do t and with 
the theoretical diffusion constant Do = /2. 

The localization length can also be determined from 
the saturation of the diffusive spreading which gives a 
localization length -^diff where ^^-^ is the time average 
of the quantum expectation value (n^) for the interval 
^max/4 < t < tmax (scc Fig. [Tj). Typically £diff is compa- 
rable to and ^in^ up to a modest numerical factor. 

In Fig. [HI we compare the three numerically calculated 
values of the localization ^in^ and £diff with the theo- 
retical expression £0 = DqT/H^ for various values of the 
classical parameters : k G {0.1, 0.15, 0.2, 0.3, 0.5, 0.7, 1.0} 
and T G {10, 15, 20, 30, 50, 70, 100} and for the fixed value 
h = 27r/(17 + 7) ^ 0.357. We see that agrees actu- 
ally very well with 2^o for a wide range of parameters and 
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FIG. 8: (Color online) Localization length determined by 
three different numerical methods as a function of the the- 
oretical value ^Q — k^T/(2h^) in a double logarithmic rep- 
resentation. The data points correspond to (red crosses) 
obtained from a linear fit of ln(|'0n|^) versus n (see Fig. [G]), 
^imp (green squares) obtained from a linear fit of (ln|'0n|^) 
versus n, -£diff (blue triangles) obtained from the satura- 
tion of quantum diffusion (see Fig. [T]). The straight full 
(blue) line represents i = 2^o- For clarity the data points 
for iimp (or ^diff) have been shifted down by a factor of 4 
(or 16). The value h = 27r/(17 + 7) 0.357 is exactly 
as in Fig. [6] while k G {0.1,0.15,0.2,0.3,0.5,0.7,1.0} and 
T G {10, 15, 20, 30, 50, 70, 100}. The Hilbert space dimension 
is mostly = 2^^ except for a few data points with the largest 
values of io where we have chosen N = 2^^. The ensemble 
averages have been performed over 100 different realizations 
of the random phase shifts. 



three orders of magnitude variation. For ii^ip and £diff the 
values are somewhat below 2£o with a modest numerical 
factor (about 1.5 or smaller) but the overall dependence 
on the parameters is still correct on all scales. 

In Fig. [HI we present a similar comparison as in Fig. [HI 
but here we have fixed the classical parameters to /c = 0.2 
and T = 50 and we vary h = 27r/(iV + 7) with 8 < TV < 
200 i. e.: 0.03132 < h < 0.7291. Again there is very 
good agreement of £^ and £in ^ with 2£q for nearly three 
orders of magnitude while for -^diff there is slight decrease 
for larger values of £q. 

The agreement £^ ^ 2io = 2DoT/h'^ in Fig. [HI for 
a very large set of classical parameters is actually too 
perfect because some of the data points fall in the regime 
where the scaling parameter x = kT^^'^ is relatively small 
(between 2.47 and 30) with a classical diffusion constant 
D well below its theoretical value I^o = k'^/2 (see Fig. [3]). 
Therefore one should expect that the localization length 
is reduced as well according to £ ^ 2DT/h^ < 2£o but 
the numerical data in Fig. [H] do not at all confirm this 
reduction of the localization length we would expect from 
a classically reduced diffusion constant. One possible ex- 
planation is that the classical mechanism of relatively 
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FIG. 9: (Color online) Same as in Fig. [8] but with fixed 
values for the classical parameters k — 0.2 and T — 50 while 
h = 27t/{N + 7) varies with the integer variable N in the 
interval 8 < iV < 200 (i. e.: 0.03132 < h < 0.7291). The 
number of realizations of the random phase shifts is 20. 

strongly correlated phases which induces the reduction 
of the diffusion constant depends on the fine structure 
of the classical dynamics of the Chirikov typical map in 
phase space, a fine structure which the quantum dynam- 
ics may not resolve if the value of h is not sufficiently low. 
Therefore the classical phase correlations are destroyed 
in the quantum simulation and we indeed observe the 
localization length i = 2io using the theoretical value of 
the diffusion constant I^o and not the reduced diffusion 
constant D. 

In order to investigate this point more thoroughly one 
should therefore vary h in order to see if it is possible 
to see this reduction of the diffusion constant also in the 
localization length provided that h is small enough. In 
Fig. [9l we have indeed data points with smaller values of 
h but here the classical parameters k = 0.2 and T = 50 
still provide a large scahng parameter kT^^'^ ^ 70.7 with 
a diffusion constant D ^ 0.9 Dq already quite close to 
Dq. In Fig. [TOl we therefore study the same values of 
h (as in Fig. [9]) but with modified classical parameters 
k = 0.2 and T = 10 such that /cT^/^ ^ 6.32 resulting 
in a diffusion constant D = 0.28591^0 well below Dq. In 
Fig. [TOl we indeed observe that the localization length 
£^ is significantly below 2£o for larger values of £0 (small 
values of h). 

We can actually refine the theoretical expression of the 
localization length in order to take into account the re- 
duction of the diffusion constant. For this, we remind 
that for short time scales t < T the initial diffusion is 
always with I^o and that only for t ^ T we observe the 
reduced diffusion constant D. We have therefore applied 
the following fit : 
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FIG. 10: (Color online) The localization length i = 
versus the theoretical value io = k'^T/{2h^) for k = 0.2, 
T = 10 and the same values of h as in Fig. [9l The 
straight full (red) upper line corresponds to ^ = 2£q. The 
lower full (blue) curve corresponds to the corrected expres- 
sion £ = 24i(4) = 24/ [/i(4) + 7/1(^0)2 + 2B^o], K£o) = 
(1 - 2A£q)/2, a ^ 0.03121, and B ^ 0.1177 due to a modi- 
fied classical dynamics (see text for explanation). The num- 
ber of random phase realizations is 50. The inset shows the 
classical diffusion for the same classical parameters k = 0.2, 
T = 10. The data points (red crosses) are (selected) nu- 
merical data and the full (blue/black) curve represents the 
numerical fit: (p^) = Do t (T + At)/{T + Bt) providing the 
two parameters A and B for the corrected localization length 
£ci in the main figure. The two straight lines correspond to 
the initial diffusion at short time scales with diffusion con- 
stant Do = k'^/2 and the final diffusion at long time scales 
with D = Do (A/B) ^ 0.2652 Do. This value of D coincides 
quite well with the value of the scaling curve in Fig. [3] at 
X = A;r^/2 6.32. 



to the classical spreading where A and B are two fit 
parameters. This expression fits actually very well the 
classical two scale diffusion with Dq for short time dif- 
fusion and with D = Do{A/B) for the long time diffu- 
sion. For k = 0.2 and T = 10 we obtain (see inset of 
Fig. ^ the values A ^ 0.03121 and B ^ 0.1177 im- 
plying D = Do {A/B) ^ 0.2652 Do which is only slightly 
below the above value D = 0.2859 Dq (obtained from 
the linear fit of the classical spreading for the interval 
lOT < t < 100 T). We can now determine a refined 
expression of the localization length using the two scale 
diffusion fit (|3T|) together with the implicit equation ([29|) 
for the critical time scale which results in the following 
equation for £ : 

.2_P2m_ i+A£ 

I -2-^-24^^^. (32) 



(31) 



This is simply a quadratic equation in £ whose positive 
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solution can be written in the form : 

24 



hi£o) + yM4F+2B4 



(33) 



with h{£o) = (1 — 2A£o)/2. One easily verifies that the 
limit D = 1^0, which corresponds to A = immediately 
reproduces ^ci(^o) = as it should be. Furthermore, the 
limit £oA > 1 (i.e. : h <C y^Do T/A) provides ^ 
£q{A/B) while for £qA <C 1 we have ^ 4 (even 

ioT A^B). 

The data points of £^ in Fig. [10] coincide very well with 
the refined expression (|33|) thus clearly confirming the in- 
fluence of the classical two scale diffusion on the value of 
the localization length as described by Eq. (|33|) . De- 
pending on the values of h the refin ed locali zation length 
£ci(^) is either given by if ^ ^ \/DqT/A or by the re- 
duced value io{A/B) if <C ^J D^T j A. In the first case 
we do not see the effect of the reduced diffusion constant 
because the value of h is too large to resolve the subtle 
fine structure of the classical dynamics. Furthermore the 
critical time scale t*, where the localization sets in, is 
below T and the momentum spreading saturates already 
in the regime of the initial short time diffusion with I^o- 
In the second case h is small enough to resolve the fine 
structure of the classical dynamics and the time scale 



is above T such that we may see the reduced diffusion 
constant D = Do{A/B) leading to the reduced localiza- 
tion length 4i(4) = £o{A/B). 



VIII. DISCUSSION 

In this work we analyzed the properties of classical and 
quantum Chirikov typical map. This map is well suited 
to describe systems with continuous chaotic flow. For 
the classical dynamics our studies established the depen- 
dence of diffusion and instability on system parameters 
being generally in agreement with the first studies pre- 
sented in [2, gI, [i3]. In the quantum case we showed 
that the chaotic diffusion is localized by quantum inter- 
ference effects giving rise to the Chirikov localization of 
quantum chaos. We demonstrated that the localization 
length is determined by the diffusion rate in agreement 
with the_general theory of Chirikov localization devel- 
oped in [6[0, H, The Chirikov typical map has more 
rich properties compared to the Chirikov standard map 
and we think that it will find interesting applications in 
future. 
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